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Abstract 

We extend a recent calculation of the nuclear energy density functional in the frame- 
work of chiral perturbation theory by computing the isovector surface and spin-orbit 
terms: {Vpp - V pnf Gd{p) + {Vpp - Vpn) ■ [Jp - Jn) Gso{p) + {Jp - Jnf Gj{p) pertaining 
to different proton and neutron densities. Our calculation treats systematically the effects 
' from Ivr-exchange, iterated Ivr-exchange, and irreducible 27r-exchange with intermediate 

CN ■ A-isobar excitations, including Pauli-blocking corrections up to three- loop order. Using 

^ \ an improved density-matrix expansion, we obtain results for the strength functions Gd{p), 

Gso{p) and Gj{p) which are considerably larger than those of phenomenological Skyrme 
forces. These (parameter- free) predictions for the strength of the isovector surface and 
spin-orbit terms as provided by the long-range pion-exchange dynamics in the nuclear 
medium should be examined in nuclear structure calculations at large neutron excess. 



PACS: 12.38.Bx, 21.30.Fe, 21.60.-n, 31.15.Ew 



^ Keywords: Nuclear energy density functional; Density-matrix expansion; Chiral pion-nucleon 
I— I ', dynamics 

> ■ 

^ 1 Introduction 

^ ' The nuclear energy density functional approach is the many-body method of choice in order 
Q to calculate the properties of medium- mass and heavy nuclei in a systematic manner [T]. Pa- 
O rameterized (non-relativistic) Skyrme functionals [21 |3l HI O [6] as well as relativistic mean-field 
models [3 [8] have been used widely and successfully for such nuclear structure calculations. 
. ^ Likewise, some constraints from chiral (pion-nucleon) dynamics and the symmetry breaking pat- 
^ \ tern of QCD at low energies have been implemented into a pertinent relativistic point-coupling 
cj ' Lagrangian in ref . [9] . 

A complementary approach in the quest for a predictive nuclear energy density functional 
[Tot fTT| [12] focusses less on the fitting of experimental data, but attempts to constrain the 
analytical form of the functional and the values of its couplings from many-body perturbation 
theory and the underlying two- and three-nucleon interactions. Switching from the conventional 
hard-core NN-potentials to low- momentum interactions is essential in this respect, because 
the nuclear many-body problem formulated in terms of the latter becomes significantly more 
perturbative. Indeed, second-order perturbative calculations provide already a good account of 
the bulk correlations in infinite nuclear matter [13j and in doubly- magic nuclei [14] . 

In many-body perturbation theory the contributions to the energy are written in terms of 
density-matrices and propagators convoluted with the finite-range interaction vertices, and are 
therefore highly non-local in both space and time. In order to make such functionals numerically 
tractable in heavy open-shell nuclei it is desirable to develop simplified approximations to these 
functionals expressed in terms of local densities and currents only. In this construction the 
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density-matrix expansion comes prominently into play as it removes the non-local character 
of the exchange (Fock) contribution to the energy by mapping it onto a generalized Skyrme 
functional with density-dependent couplings. For some time the prototype for that has been 
the density- matrix expansion of Negele and Vautherin |15j, but recently Gebremariam, Duguet 
and Bogner [16] have developed an improved version for spin-unsaturated nuclei. They have 
demonstrated that phase-space averaging techniques allow for a consistent expansion of both the 
spin-independent (scalar) part as well as the spin-dependent (vector) part of the density-matrix. 
The accuracy of the new phase-space averaged density-matrix expansion and the original one 
of Negele and Vautherin [15j has been gauged via the Fock energy (densities) arising from 
(schematic finite-range) central, tensor and spin-orbit interactions for a large set of semi-magic 
nuclei. For a central force the Fock energy depends primarily on the spin-independent (scalar) 
part of the density-matrix and a few percent accuracy is reached for both variants of the density- 
matrix expansion. On the other hand the Fock energy due to a tensor force is determined by 
the spin- dependent (vector) part of the density-matrix. In that case the original density-matrix 
expansion of Negele and Vautherin [TJ] leads to an error of about 50%, whereas the new one 
based on phase-space averaging techniques reduces the error drastically to only a few percent. 
For further details on these extensive and instructive test studies we refer to ref. [16j. 

In order to match with these new developments the nuclear energy density functional as 
it emerges from chiral pion-nucleon dynamics has been recalculated recently in ref. [17J. That 
calculation has treated for isospin-symmetric (i.e. N = Z) nuclear systems the effects from 
Ivr-exchange, iterated Ivr-exchange, and irreducible 27r-exchange with intermediate A-isobar 
excitations, including Pauli-blocking corrections up to three-loop order. It has been found that 
the effective nucleon mass M*{p) entering the energy density functional becomes identical to 
the one of Fermi-liquid theory when employing the improved density-matrix expansion. The 
strength F^{p) of the (Vp)^ surface term as provided by the pion-exchange dynamics comes 
out in good agreement with empirical values in the density region po/2 < p < po = 0.16 fm"'^. 
The spin-orbit coupling strength -Fso(p) receives contributions from iterated Ivr-exchange (of 
the "wrong sign") and from three- nucleon interactions mediated by 27r-exchange with virtual 
A-excitation (of the "correct sign"). In the region around po/2 — 0.08 fm~^ where the spin- 
orbit interaction in nuclei gains most of its weight these two components tend to cancel, thus 
leaving all room for the short-range spin-orbit interaction. As a matter of fact the empirical 
value Fi^'^P) ~ 90 MeVfm^ of the spin-orbit coupling strength in nuclei agrees well with the 
one extracted from the short-range spin-orbit component of any realistic NN-potential [18]. 
This part of the NN-interaction drives at the same time the strong Lorentz scalar and vector 
mean-fields [19J on which the whole success of the relativistic Dirac phenomenology rests. 

The purpose of the present paper is the extend the calculation of the energy density func- 
tional in ref. [IT] to isospin-asymmetric many- nucleon systems with different proton and neutron 
densities. The supplementary isovector surface and spin-orbit terms play an important role in 
the description of long chains of stable isotopes and for nuclei far from stability. Our paper 
is organized as follows. In section 2, we recall the improved density matrix expansion of Ge- 
bremariam, Duguet and Bogner [TB] whose Fourier transform to momentum space provides 
the adequate technical tool to calculate the nuclear energy density functional in a diagram- 
matic framework. In section 3, we present the analytical results for the density-dependent 
strength functions Gd{p), Gso{p) and Gj{p) belonging to the isovector surface and spin-orbit 
terms. These analytical expressions give individually the effects due to Ivr-exchange, iterated 
Ivr-exchange, and irreducible 2vr-exchange with intermediate A-isobar excitations, including 
Pauli-blocking corrections up to three-loop order. Section 4 is devoted to a discussion of our 
numerical results and it ends with some concluding remarks and an outlook. 
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2 Energy density functional and improved density-matrix 
expansion 



Let us begin with writing down the exphcit form of the isovector surface and spin-orbit terms 
in the nuclear energy density functional: 

^iv[Pp,Pn, ^p, ^] = Pp-^ Pnf Gd{p) + {V Pp-S/ Pn)-{Jp- Jn)Gso{p) + {Jp- Jnf Gj{p) . (1) 

Here, 

Pp,n{r) = %P = E ^i:i\n^^;:lir) , (2) 

are the local proton and neutron densities written in terms of the corresponding (local) proton 
and neutron Fermi momenta kp^nir), and expressed as sums over the occupied single-particle 
orbitals The spin-orbit densities of protons and neutrons are defined similarly: 

JpAf) = E ^Ifin^ X vvi/(5(f ) . (3) 

a 

Furthermore, Gd{p), Gso{p) and Gj{p) in eq.(l) denote the associated strength functions de- 
pending on the total nucleon density p = Pp+Pn- In Skyrme parameterizations [2l[3llll[5l|6] these 
are just constants, Gf^'^ = -[3ti(2xi + 1) + ^2(2x3 + l)]/64, GfJ"^ = Wq/A, Gf^^ = (ti -t2)/32, 
whereas in our calculation their explicit density-dependence originates from the finite-range 
character of the Ivr- and 27r-exchange interaction. 

The starting point for the construction of an explicit nuclear energy density functional 
£^iv[pp, Pn, Jp, Jn] is the bilocal density-matrix as given by a sum over the occupied energy eigen- 
functions: J2a ^p"2(^~^/2)^p°2^(^^+^i/2). According to Gebremariam, Duguet and Bogner [TB] 
it can be expanded in relative and center-of-mass coordinates, a and r, with expansion coeffi- 
cients determined by local quantities (nucleon density, kinetic energy density and spin-orbit den- 
sity). As outlined in section 2 of ref.[l7j the Fourier transform of the expanded density- matrix 
with respect to both coordinates a and f defines in momentum space a "medium- insertion" 
r{p,q) for the inhomogeneous many- nucleon system. It is straightforward to generalize this 
construction [17J to the isospin-asymmetric situation with different proton and neutron local 
densities Pp,n(^) and Jp,„(r). We display here only that part of the medium-insertion T{p,q) 
which is actually relevant for the diagrammatic calculation of the isovector surface and spin- 
orbit terms introduced in eq.(l): 

-j^Sikf-\p\)T3a-\px (Jp- J„)] + ... j. (4) 

When working to quadratic order in deviations from isospin symmetry (proton-neutron differ- 
ences) it is sufficient to use an average Fermi momentum kf in the prefactor of the isovector 
spin-orbit density Jp — Jn- The double line in the left picture of Fig. 1 symbolizes this medium 
insertion together with the assignment of the out- and in-going nucleon momenta p±q/2. The 
momentum transfer q is provided by the Fourier components of the inhomogeneous (matter) 
distributions Pp,n{r) and Jp,n{r). Note that in comparison to the version of T{p,q) which fol- 
lowed from Negele and Vautherin's density-matrix expansion [15] the weight function of the 
spin-orbit densities Jp^nif) has changed from 5{kf — \p\) — kfS\kf — \p\) to —36{kf — \p\). 
Furthermore, a pairwise filling of time-reversed orbitals a for both protons and neutrons has 
been assumed, so that (various possible) time-reversal-odd fields do not come into play [1]. 
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Figure 1: Left: The double-line symbolizes the medium insertion defined by eq.(4). Next are 
shown: The one-pion exchange Fock diagram and the iterated one-pion exchange Hartree and 
Fock diagrams. 

3 Diagrammatic calculation 

In this section we present analytical formulas for the three density- dependent strength functions 
Gd{p), Gso{p) and Gj{p) as derived (via the improved density-matrix expansion pJlj) from 
Ivr-exchange, iterated Ivr-exchange, and irreducible 27r-exchange diagrams with intermediate 
A-isobar excitations, including Pauli-blocking corrections up to three-loop order. We give for 
each diagram only the final result omitting all technical details related to extensive algebraic 
manipulations and solving elementary integrals. In essence the calculation of ref. [IT] gets just 
modified by relative isospin factors —1/3, ±2/3, —5/3 occurring at various places. 

3.1 One-pion exchange Fock diagram 

The non-vanishing contribution from the Ivr-exchange Fock diagram shown in Fig. 1 reads: 



gAp) 



(32m,/,)2u6 



Au^ - - ln(l + Au^) 



(5) 



where we have introduced the convenient dimensionless variable u = kf/m.,^. This contribution 
to Gj{p) is just —1/3 of the contribution to the isoscalar strength function Fj{p) (see eq.(ll) 
in ref. ^j) as a consequence of the isospin trace: tr[ra(Js + T3,J^)Ta{Js + T'iJv)] = QJs ~ 2J^. 



3.2 Iterated l7r-exchange Hartree diagram with two medium inser- 
tions 

The two-body contributions from the iterated Ivr-exchange Hartree diagram in Fig. 1 read: 

Gsoip) = ^t^^Iqu^ + I ln(l + 4^2) - 8m arctan 2u] , (7) 

which are —2/3 of the respective isoscalar contributions [T71 |2U]. Let us briefiy explain the 
mechanism which generates the strength function Gd{p)- The exchanged pion-pair transfers 
the momentum q between the left and the right nucleon ring. This momentum q enters both 
the pseudovector vrA^-interaction vertices and the pion propagators. After expanding the inner 
loop integral to order the Fourier transformation in eq.(4) converts this factor into a 
factor (Vfcp - Vknf ~ 7r^(Vpp - Vpnf/k). The rest is a solvable integral over the product of 



4 



two Fermi spheres of radius kf. The isovector spin-orbit strength function Gso{p) arises from 
the spin-trace: tr[a ■ (/ + q/2) a ■ {I — q/2) a ■ {pi2 x J^)] = 2i {qx I) ■ (pi 2 x Jv), where q gets 
again converted to V/cp — V/cn — 7r^(Vpp — Vpn)/k'j by Fourier transformation. 



3.3 Iterated l7r-exchange Fock diagram with two medium insertions 

We find the following contributions from the right Fock diagram in Fig. 1 with two medium 
insertions on non-neighboring nucleon propagators: 

5g\M r4_ , ^ . . 1 , (1 + ^2^2 



^ , , 5g\M (A. 

Gd(p) = 7 — — -7< — (2 arctan2n — arctanw) 

37rm^(8/^)* [u 



■\ In 



;1 + 2m2)(i + 4m2) 



4 



1 + 2^2 Jo 



+ — I dx 



3 + 18x^ + 16x^ r 



arctan x — arctan 2x 



(8) 



hg\M 



2TTm^{Af^u) 



+ dx 

Jo 1 + 2x2 



X 



^1 -|- 4a; ) arctan 2a: — 4a; (1 + x ) arctan x 



(9) 



15g\M [2 2 /■« 



1 + 2x2 

+2(6x^ - 2u^x'^ + 6x^ - 2u^x^ + 2x^ - m^) arctan x 



■u^)(l + 4x^) arctan 2x 



(10) 



One notices a relative isospin factor of —5/3 in comparison to the respective isoscalar contribu- 
tions [171 120] which comes from the isospin trace tr[ra( + T3Jv)TbTa{Js + 'T3Jv)'''b] = 10 J2 — 6 J2 
of this diagram. 



3.4 Iterated Irr-exchange Hartree diagram with three medium in- 
sertions 

In our way of organizing the many-body calculation, the Pauli-blocking corrections are repre- 
sented by diagrams with three medium insertions. The corresponding contributions from the 
iterated Ivr-exchange Hartree diagram shown in Fig. 1 read: 

2g\M fl f\ , l + yr4M2|/2(6 + 51^2^2 + 92mV) 
dy m 



G,{p) 



37r2m7r(4/7r)'' [u Jo 
-21n(l + 4MV 



+ I dx — 

J-l 



l-yl 3(1 + 4^2^2)3 

x2 /•! , s^s'(3 + 9^2 + As^) , u + xy 
dy — In 



l + s 



2\4: 



u — xy 



(11) 



G,o{p) 



2g\M 



X 



TT^m^iAUyJo '^''u^J-fi 



, u + xy M(5x2y2 _ ^u"^) 



— x2|/2 



X 



+ 



5s + 



'l + s 



2\2 



6 arctan s 



2sh' 



Gj{p) 



(1 + S2)2 

?,g\Mu^ 



4m — (s + 2xy) In 



u — xy 

us^iu^ + x2?/2) 

1 + ^2)2(^2 — a;2y2) 

U + xy' 



u — xy 



dy 



Wm^f^ Jo ^ (1 + 4^2^2)2 



2y + {y^-l)\n 



1 + y 



[12) 



(13) 



with the auxiliary function s = xy + y/u'^ — x'^ + x'^y'^ and its partial derivative s' = uds/du. 
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3.5 Iterated l7r-exchange Fock diagram with three medium inser- 
tions 



The evaluation of this diagram is most tedious. It is advisable to split the contributions to the 
strength functions Gso(p) and Gj{p) into " factorizable" and "non-factorizable" parts. 

These two pieces are distinguished by the feature of whether the nucleoli propagator in the 
denominator can be canceled or not by terms from the product of vrA'^-interaction vertices in 
the numerator. We find the following "factorizable" contributions: 



Gd{p) 



9aM 



95 arctan 2u 1 + 6m 



67r2m, (4/^)4 [ 8(1 + ^2) 
9 + IOOm^ + 184^4 - 192^6 



+ 

+ ^ / dx 
Jo 

X — 4u 



32u^ 
ln(l + 4m^ 



7u 



1 + 4^2 



x^ 



x^ 



X -\- Au 



+ 



2x 



_ ]_ 

4m+ — (l + M^) 
x^ 

2x 



l+{u + x)^ l + {u-x)^ [l + (M + a;)2]2 [1 + (m - a;)2]2 



(14) 



Gsoip) = 



9\M 



3 + 14^2 + 2m^ 



19m^ - ln(l + 4m^) arctan 2m - 60m^ - 8m 

J M 



+ 



+20 



2, 3 + 20m2 + 16m\ 2/. . 2, 
2.3 ) '"'^ + *" 



ln(l + 4m2) 



lo'^i 



M^ - 2m' - 3 - 3x-'(l + m')'' + (3 + 7m')x' - 5x 



+6ux-^{l + u^fL - 3u^x-^{l + M^) 1 1 , 



(15) 



Gj{p) = 



3g\M { 185m _ 

7r2m^(8/,M)4^~ 
1 99m2 - 8 - 68m^ 



H h 

M 



16m3 



4m — 35 arctan 2m + 



ln(l + 4m^) + 5 dx 



1 + 4m2 
16m5 



L2 



M^ 



In2(l + 4M2) 



2\4 



2x2 



(l + M^) 



+2(m^ - 1)(1 + M ) - (5 + 2m' + 5M^)a;' + (6 + 10M')a;^ - 



llx^ 



+- 



M 



1 + M^)^ - 3m^ - 2m^ + 1 



x^ 



2\2 I 



2a;2 



(16) 



with the logarithmic function: 



1 l + (M + a;)2 

L{x, m) = — In —— ^ 

Ax l + [u- xy 



(17) 



The "non-factorizable" contributions (stemming from nine-dimensional principal value integrals 
over the product of three Fermi-spheres of radius kj) read on the other hand: 



Gdip) 



9\M 



\ f\ yz9{y^ + -1) (2uy^9{y)9{z) 



-2mV(3j_4^^ 
. (1 + 4m2^2)2 



dy / dz - — |— p= 
1 J-i \yz\y/y 



ln(l + Au^z^) 



2 + ,^2 _ 1 \ 1 + 4^2^^2 

« bx^ss'tH'il ~ 5s2 - 2s^) 



-\- dr 

Jo 4m6(1 + s2)3(1 + ^2) 



, (18) 



Gso{p) 



9aM 



X 2uz — arctan(2M2;) 



/ I — 1/2,2 f 

1 J-i \yz\y/y^ + — 1 



/o 2u8(l + s2)(l+t2 



ly^z9{y)9{z) 

1 + 4M2y2 

{st + 5X2; — te?/) \ , (19) 



dy 



z 



4:U^z^ - ln(l + Au^z^ 



-1 \yz\y/]f~+^ 



X 



1 U(l +4M2y2) 
5x'^s3s't3t'(?/2 + ^2 _ 



-I- / dr 

/o 4MlO(l + s2)(l+t2) 



(20) 



with the auxihary function t = xz + y/u'^ — x'^ + x'^z'^ and its partial derivative t' = udt/du. 
For the numerical evaluation of the d2;-double integrals in eqs. (18, 19,20) it is advantageous 
to first antisymmetrize the integrands in y and z and then to substitute z = y/y'^C^ + 1 — y'^. 
This way the integration region becomes equal to the unit-square < C < 1. 



3.6 Irreducible two-pion exchange 

At next order in the small momentum expansion comes the irreducible 27r-exchange including 
(also) intermediate A-isobar excitations. We employ a (subtracted) spectral-function represen- 
tation of the TrA^A-loops and find the following non-vanishing (two-body) contribution: 



gAp) 



— / rf/i Im(Vc - Wc) 



+Im(VT - Wt) 



3/i ^ 2k^ 



f 



P^ / 2 



Akj) In ( 1 + 



P' 



4 



(21) 



The imaginary parts ImV^, ImPFc", ImV^ and ImPVV of the isoscalar and isovector central and 
tensor NN-amplitudes due to 27r-exchange with single and double A-excitation can be found in 
section 3 of ref.[2T]. The additional contributions from the irreducible 27r-exchange with only 
nucleon intermediate states are accounted for by inserting into eq.(21) the imaginary parts: 



ImW, 



c 



P' 



Ami 



3vr/i(4^) 



4m2 (1 + Ag\ - 5g\) + fi\23g\ - lOgl - 1) + 



A8g\ml 



4m2 



ImVr 



(22) 



(23) 



vr/i(4/.)4 

At leading order the irreducible 27r-exchange generates no spin-orbit NN-interaction. It emerges 
first as a relativistic 1/M-correction. In order to see the size of such relativistic effects we have 
evaluated the energy density functional with a two-body interaction composed of the (isoscalar 
and isovector) spin-orbit NN-amplitudes Vso and Wso written in eqs. (22,23) of ref.[22]. We 
find with it the following contribution to the isovector spin-orbit coupling strength: 



Gsoip) 



7rM(4/^)4 I 15u^ 



^ "^^ln(l + «^)-n^]+^(ll,i-4) 



+ 



4 - 10g\ 4m , 

6U 5 



arctan u 



(24) 



which has been subtracted at p = in order to eliminate (regularization dependent) short- 
distance components. As a consequence of that subtraction only the Fock terms are included 
in the expressions in eqs. (21, 24). 
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Figure 2: Hartree and Fock three-body diagrams related to 27r-exchange with virtual A-isobar 
excitation ("Fujita-Miyazawa [23] mechanism"). 



3.7 Three-body diagrams with A-excitation 

The Pauli-blocking correction to the 27r-exchange with single A-excitation is equivalent to the 
contribution of a (genuine) three-nucleon force. In fact, one is dealing here with the same 
three-nucleon interaction as originally introduced by Fujita and Miyazawa [23]. The pertinent 
Hartree and Fock three-body diagrams related to 27r-exchange with virtual A-excitation are 
shown in Fig. 2. Returning to the medium insertion written in eq.(4) we find from the left 
three-body Hartree diagram in Fig. 2 the following (non- vanishing) contribution: 



Gj{p) 



TV 



=A(4^) 



6 . 
4m - 

u 



6u 



l + 4u2 



-^Ml + 4..=)} 



(25) 



with A = 293 MeV the delta- nucleon mass splitting. We have used the value 3/a/2 for the ratio 
between the ttA^A- and vrA^A^-coupling constants. The vanishing of contributions to Gd{p) and 
Gso{p) has the following reason. The pertinent isospin trace generates an expression that is odd 
under exchange of the momenta pi and p2 attached to the (upward and downward running) 
nucleon lines of the (left) closed nucleon ring. The remaining factor from the two-pion exchange 
interaction is however even (under pi ■v^ P2) and so the whole expression integrates to zero. 

— * — * 

The three-body effects on the energy density functional £^iv[Pp,Pn, Jp, Jn\ are completed by 
the contributions from the central and right Fock diagrams in Fig. 2, which read: 



Gd{p) 



112m 



37r2A(8/^)4 [ 1 + 4m2 

64m 14 20 

3 



16^ 



20 ^ 30 



u 



+ 



15 



1 + 4m2 



u 



Vn(l 



+ 4m^ 



r(3 + IOm^ + Sm^) In2(l + 4M^ 



(26) 



Gso{p) 




5 
4^ 



- — \ ln(l + 4u' 
3ti/ 



— ^—mu^ + 80m^ + 36^2 + 5) \^{\ + 4^2) I 
64^^^ J 



(27) 



gAp) 




ln(l + 4m^ 



16 )'"(^+^"''-4.3 



arctan 2u 
3 1 



u 



8 



787m 3935^3 1696^5 3 + IGu^ + 144m^ , „^ „^ 

H 11^(1 + Au^) 

16 48 15 64m7 ^ ^ ^ 



+ dx 



" . r 3L 



L 



^(1 + n^)^(3..2 - 1) - g(l + «^)^ + (1 + n^)^ 



X (94^2 - 145 - 257u^) + 4x^(117m^ + 35u^ - 145m^ - 63) - 103a;' 
+x\lS2u'^ - 165 - 597m^) + 2x^(199^2 - 77)^ ^ 



4m 



105,^ 2N5 



7 



+ — (1 + uy{3 - 49u^) + 18(1 + u^){37u^ - + 23) 



105 



i + «T + ^(i + «y(3 + ii«^)}} 



(2J 



with L(x, u) defined in eq.(17). A good check of all formulas collected in this section is provided 
by their Taylor series expansion in kf. Despite the superficial opposite appearance the leading 
term in the ^/-expansion is always a non- negative power of kf (which is higher for three-body 
contributions than for two-body contributions). 

4 Results and discussion 

In this section we present and discuss our numerical results obtained by summing the series of 
contributions given in section 3. The physical input parameters are: qa = 1-3 (nucleon axial 
vector coupling constant), = 92.4 MeV (pion decay constant), = 135 MeV (neutral pion 
mass) and M = 939 MeV (nucleon mass). We recall that with these physical parameters and 
a few adjustable short-distance couplings the nuclear matter equation of state E{p) and many 
other nuclear matter properties [21] can be well described by the chiral pion-nucleon dynamics 
treated up to three-loop order. 

Fig. 3 shows the strength function Gd{p) belonging to the isovector surface term (Vpp— Vp„)^ 
plotted versus the nucleon density p = 2fcj/37r^. One observes that the leading result due to 
the iterated Ivr-exchange (eqs. (6,8,11,14,18) shown separately by the dashed line) is almost 
not changed by the inclusion of three-body contribution eq.(26) related to 27r-exchange with 
virtual A-excitation. Moreover, the density dependence of Gd{p) is rather weak in the entire 
region < p < 0.2 fm~^. The horizontal lines in Fig. 3 correspond to the (constant) values 
Cf^^ = — [3ti(2xi + 1) + t2{2x2 + l)]/64 of three phenomenological Skyrme forces SkM* [5], 
SkP [Ij and Sly [5]. These values are of the same negative sign, but considerably smaller 
in magnitude than our parameter-free prediction resulting from the long-range pion-exchange 
dynamics in the nuclear medium. It remains to be seen how well these (predicted) much larger 
negative values of Gd{p), which energetically favor large differences in the density-gradients of 
protons and neutrons, perform in actual nuclear structure calculations at large neutron excess. 

Of particular interest is the strength function Gso[p) of the isovector spin-orbit coupling term 
(Vpp — Vpn) ■ ( Jp — Jn) as provided by the explicit pion-exchange dynamics. As one can see from 
Fig. 4 the leading result due to the iterated Ivr-exchange (eqs. (7,9,12,15,19) shown separately by 
the dashed line) is again not changed by the inclusion of the 27r-exchange and associated three- 
body contributions eqs. (24,27). This feature is markedly different from the isoscalar spin-orbit 
coupling strength -Fso(p) where an almost complete cancellation between these two components 
has occurred around po/2 = 0.08 fm^^ (see Fig. 5 in ref.[T7j). The basic reason for this totally 
different behavior is the absence of a strong isovector three-body spin-orbit coupling arising 
through the Fujita-Miyazawa mechanism (i.e. from the dominant three-body Hartree diagram 
in Fig. 2). There is of course in addition the short-range spin-orbit interaction which has to 
account practically for the full strength -Fj^^P^ ~ 90 MeVfm^ in the isoscalar channel and which 
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Figure 3: Strength function Gd{p) of the isovector surface term (Vpp — V PnY in the nuclear 
energy density functional versus the nucleon density p = 2A;j/37r^. Dashed line: Iterated Itt- 
exchange only. Full line: 27r-exchange and associated three-body contributions added. 
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Figure 4: Strength function Gso{p) of the isovector spin-orbit couphng term (Vpp — Vp„) • {Jp — 
Jn) in the nuclear energy density functional versus the nucleon density p = 2k^/3n'^. Dashed 
line: Iterated l7r-exchange only. Full line: 27r-exchange and three-body contributions added. 
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Figure 5: Strength function Gj{p) belonging to the squared isovector spin-orbit density (Jp — 
JnY in the nuclear energy density functional versus the nucleon density p = 2A;j/37r^. Dashed 
line: Ivr- and iterated Ivr-exchange only. Full line: 27r-exchange and three-body contributions 
added. 

contributes also in the isovector channel with a reduced weight 1/3. Combined with our result 
for Gso{p) from the long-range pion-exchange dynamics one would then have a situation where 
the isoscalar and isovector spin-orbit coupling strengths are about equally strong. As discussed 
in ref.[25] the isotope shifts of the charge radii in the Pb region (i.e. around the shell-closure 
N = 126) can provide a sensitive test for strength of the isovector spin-orbit coupling. However, 
according to their limited analysis (in the vicinity of the stability line) no definite choice could 
be made between several density- dependences of the neutron spin-orbit potential: ~ Pp + 2pn 
for the Skyrme force, ~ Pp + Pn for the relativistic mean-field model, or ~ pp for the generalized 
functional SkI4 [25]. It remains to be seen whether the proportionality of the neutron spin-orbit 
potential to the gradient of the neutron density Vpn as suggested by the present calculation 
leads to results which are in accordance with the existing precise experimental data. 

Finally, we show in Fig. 5 the strength function Gj{p) belonging the squared isovector 
spin-orbit density (Jp — J„)^ in the nuclear energy density functional as a function of the 
nuclear density p. One observes that the inclusion of the subleading 2'7r-exchange (two- and 
three-body) contributions eqs. (21,25,28) significantly affects the outcome for Gj{p). For better 
orientation, we reproduce also by the dashed-dotted line in Fig. 5 the leading contribution 
from the Ivr-exchange Fock diagram alone (see eq.(5)). Interestingly, in the density region 
around po/2 = 0.08 fm~^ the Ivr-exchange approximation and the complete result (full line) 
agree roughly with each other. The constant values G^j^^ = (ti — 12)/32 from phenomenological 
Skyrme forces [31 HI [5] are of similar magnitude but opposite in sign. Note that the strength 
function Gj{p) comprises in particular the non-local Fock contributions from tensor forces, 
whose long-range isovector component is determined by the Ivr-exchange. This correspondence 
gives a rough understanding of the features visible in Fig. 5. Another interesting side effect of 
the (Jp — Jp)^ term in the nuclear energy density functional is that it gives rise to an extra 
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spin-orbit mean-field 2Gj{p) {Jp — J„) in addition to the "normal" one, Gso{p) {S Pp — ^Pn)- 
It would be interesting to investigate the role of this additional (nucleus-dependent) spin-orbit 
mean- field together with the density-dependence of Gj{p) as predicted by in medium chiral 
perturbation theory. 

In summary, we have used the improved density- matrix expansion of ref . [16] to calculate the 
strength functions of the isovector surface and spin-orbit terms in the nuclear energy density 
functional as provided by the long-range pion-exchange dynamics in the nuclear medium. These 
predictions together with the ones in ref. [1^ for the isoscalar strength functions should be 
examined and explored in nuclear structure calculations at small and large neutron excess. 
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